The effect of hidden color channels on Nucleon-Nucleon interaction* 
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This letter reports the nucleon-nucleon(A^7V) interaction obtained from multi-channel, including 
hidden color channels, coupling quark model calculation. The results show that the hidden color 
channels coupling provides the intermediate range attraction which is usually assumed to be due 
to multi-TT or a meson exchange and that the short and intermediate range A^'A^ interaction can be 
described solely by the fundamental quark-gluon degree of freedom of QCD. 
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The TV TV interaction has been studied for more than 70 
years since the discovery of neutron and the meson ex- 
change model proposed by H. Yukawa. The effective one 
boson exchange model (OBE) fT, ^ and the chiral per- 
turbation theory (ChPT) [3] fit the NN interaction data 
below the tt production threshold quantitatively well but 
with ~30 (ChPT) to more than 40 (OBE) adjustable 
parameters. Quark models [1, Q fit the NN interac- 
tion data but quantitatively not as well as OBE and 
ChPT. However they use much less adjustable param- 
eters, 9 for [5, and fit the hadron spectroscopy simulta- 
neously. Lattice QCD has got the static qq interaction 
of 2, 3, 4 and 5 quark systems and the qualitative NN 
interaction in the quenched approximation recently Q. 
The static interaction can be expressed as color Coulomb 
plus the linear confinement proportional to the minimum 
length of the color fiux with junctions. Unquenched lat- 
tice results only modify the linear confinement with color 
screening but have not found the trace of meson exchange 
yet 0. One expects lattice QCD will describe NN inter- 
action by means of the fundamental quark-gluon degree 
of freedom of QCD directly. 

The short and even the intermediate range NN inter- 
action should be related to the nucleon internal structure 
as emphasized by P.W. Anderson Q. The long standing 
fact that the "^Si and -^Sq NN interactions are similar 
to the hydrogen molecular forces in the spin singlet and 
triplet states respectively except the energy and length 
scale difference, which certainly calls for an explanation 
but has never been explained in OBE and ChPT. Quark 
model should be able to do this job, because the NN in- 
teraction is treated as a remnant of the fundamental color 
interaction due to gluon exchange between color quarks 
of two color singlet nucleons, which is quite the same as 
the hydrogen molecular force is a remnant of the fun- 
damental electro-magnetic interaction due to photon ex- 
change between charged particles of two electric neutral 
hydrogen atoms. However it is well known that the effec- 
tive one gluon exchange Breit-Fermi interaction plus two 



body confinement quark models only obtained the NN 
short range repulsion but can not describe the interme- 
diate and long range attraction. Chiral quark model [3| 
introduced the tt and a meson exchange to provide the 
intermediate and long range attraction. We proposed a 
model 0] in which we introduced the quark mutual delo- 
calization (or percolation) to describe the spacial mutual 
distortion of interacting nucleons and a phenomcnologi- 
cal color screening to model the color space mutual dis- 
tortion (hidden color channel effect). In this report we 
show the NN scattering results obtained from the multi 
channel coupling calculation with hidden color channels 
coupling directly. These results show that the hidden 
color channels coupling does give rise to the NN inter- 
mediate range attraction and that the short and inter- 
mediate range NN interaction can be described by the 
fundamental quark-gluon degree of freedom directly. The 
mechanism of the short and intermediate range NN in- 
teraction is quite the same as that of the hydrogen molec- 
ular force and so provide a natural explanation of the 
similarity between nuclear and molecular forces. 

The NN system (a six valence quark system) can be 
consisted of two color singlet nucleons as in the usual 
hadron degree of freedom description, but also of two 
color octet nucleons coupled to an overall color singlet 
six quark state as shown in Fig. 1. The latter is called 
hidden color channel and because of color confinement, 
these hidden color channels exist in the two nucleon over- 
lap region only. The quark cluster model channel wave 
function can be expressed as. 



^{NN) = AmNi)ij{N2)]cisF{R)], 



(1) 
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where A is the anti-symmetrization operator, ip{Ni) 
is the nucleon internal wave function, Xct is the color 
part which can be color singlet or octet, CiiSi is the 
SU{4) D SU^{2) X 3^(2) spin-isospin part, V(l») is 
the orbital part and is the internal Jacobian coor- 
dinates. In order to simplify numerical calculation we 
assume it is a production of Gaussian function with a 
size parameter b. [■ ■ ■]cis means coupling the individ- 
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ual nucleoli color-isospin-spin into overall color singlet, 
total isospin / and spin S with SU^S), SW {2), SW^ (2) 
Clebsch-Gordan Coefficients. F{IV) is the relative orbital 
wave function and R is the relative coordinate between 
two nucleon center of mass coordinates, R — — R-Af2- 
The quark delocalization description can be done as the 
same in [sj and we will not repeat here to save the space. 



Table 1. Quark model parameters. 




Fig. 1. color structures of NN system. 

The model Hamiltonian of the NN system is chosen 
to be. 
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Here Sij is quark tensor operator, Y{x) and H{x) are 
standard Yukawa functions, Tc is the kinetic energy of 
the center of mass. All other symbols have their usual 
meaning. 

The model parameters are fixed as follows: the u, d- 
quark mass difference is neglected and mu=md is as- 
sumed to be exactly 1/3 of the nucleon mass, namely 
TO.,(=m(i=313 MeV, the tt mass takes the experimental 
value, the A^r takes the same values as in Ref.[4|], namely 
A7r=4.2 fm^^. The chiral coupling constant ach is deter- 
mined from the ttNN coupling constant as usual. The 
rest parameters Oc, Vq, and as are determined by fitting 
the nucleon and A masses and the stability of nucleon 
mass with respect to the variation of the nucleon size b. 
All parameters used are listed in Table 1. 



quark masses 


(MeV) 


313 


nucleon size 


b (fm) 


0.518 




TOtt (fm"-'^) 


0.71 


TT 


A^ (fm-i) 


4.2 




ach 


0.027 


Confinement 


adMeV fm-2) 


56.75 




Vo (fm2) 


-1.3590 


OGE 


as 


0.485 



This study is focused on the hidden color channel ef- 
fect and so in order to simplify the numerical calculation 
we approximate the nucleon internal orbital wave func- 
tion by a product of two Gaussian functions of the Ja- 
cobian coordinates. This approximation will make the 
exchange interaction matrix elements decreasing with R 
as /"^^ and so definitely can not describe the one tt 
exchange Yukawa interaction. Therefore we include the 
one TT exchange in our model Hamiltonian. The many 
body linear confinement interaction as obtained in lattice 
QCD is also approximated by two body quadratic one as 
shown in Eq.(2) and the color screening due to qq exci- 
tation observed in unqucnchcd lattice QCD is neglected 
too. These approximations will hinder the quark mutual 
delocalization and quite possible reduce the effective NN 
attraction as otherwise due to quark delocalization. All 
of these shortcomings are left for further improvement. 

The Kamimura variational method |8] is used to do 
multi channel scattering calculation. The channels in- 
cluded in different partial waves are listed in Table 2. 
The calculated partial wave phase shifts and the SP07 
data points are shown in Fig. 2-5. The direct extension 
of the color dependent two body confinement model can 
not describe the NN scattering quantitatively well even 
after including hidden color channels coupling as can be 
seen from the first set results (the solid lines in Figs. 2-5). 
The color dependent two body confinement interaction 
is consistent with the lattice QCD results only for two 
and three quark systems in color singlet states but in- 
consistent with the many body interaction obtained in 
lattice QCD ones for multi-quark systems For multi- 
quark systems and color octet nucleons, quark pairs are 
not always in color antisymmetric state but also color 
symmetric ones. The color factor A^- \j will give rise to 
anti-confinement interaction for symmetric quark pairs 
[To| . There is no sound theoretical reason to extend 
the color dependent two body confinement interaction 
to multi-quark system. We adjust the confinement in- 
teraction strength Oc for every partial waves to fit the 
NN phase shifts of SP07 by two recipes, recipe 1: ad- 
justing the flc to fcflc for the coupling between hidden 
color channels and the color singlet channels only; recipe 
2: adjusting the confinement strength Cc to kac not only 
for the coupling between hidden color channels and the 
color singlet ones but also for the hidden color channel 
themselves. The best fit of the multiplicative factor fc's 
are listed in Table 2 too. One can see from the dotted 
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(recipe 1) and dashed (recipe 2) lines in Figs. 2-5 that 
by including the hidden color channels and adjusting the 
color confinenient interaction strength Ac, both adjusting 
recipes can fit the NN scattering phase shifts quantita- 
tively. We take these results as an indication that the 
short and intermediate range NN interaction can be de- 
scribed solely by the fundamental quark-gluon degree of 



freedom. It is the nucleon internal structure and its dis- 
tortion both in orbital and color spaces which give rise 
to the NN short range repulsion and intermediate range 
attraction and these are quite the same as the atomic in- 
ternal structure and its distortion in orbital and electric 
charge spaces which give rise to the hydrogen molecular 
interaction. 



Table 2. The channels used in NN scattering calculations and the confinement strength for each channel. 



T J 
1 J 


channels recipe k 


1 


^SiCDi) : NN, AA, ^Ag ^As,-^ Ns ^Ns, ^iVg ^Ns, ^Ns ^Ns 1 1.40 
L>i : A/1, iVg iVg z i.6o 


iPi : A^iV, AA, 2Ag ^Ag, ^Ns ^iVg, ^Afg ^Ns 1 1.80 


2 


: NN, AA, ^Ag ^Ag, ^Ns ^iVg, ^Ns ^iVg, ^Ns ^TVg 1 1.00 

U2 ' AiA, -iVg -/Vg Z i.UU 


3 


3Z?3 : NN, AA, ^Ag ^Ag, ^N^ ^iVg, ^iVg ^N^, ^ATg ^Ns 1 2.40 
"^SsCDa) : AA, ^Ns ^Ns 2 2.20 


1 


1^0 : A^A^, AA, ^Ag ^Ag, ^Ns ^Ns, ^iVg ^Ag, ^TVg ^Ns 1 1.42 
^Do : ATA, AA, ^Ns ^Ag, ^Ng ^Ng, ^A^g ^ATg 2 1.39 


: NN, NA, AA, ^Ag ^Ag, ^Ns ^Ag, ^Ng ^Ng, ^Ng ^Ns, 1 1.10 
^iVg 2Ag, ^Ns ^Ns 2 1.10 


1 1 


'■'Pi : NN, N A, AA. ^Ag ^Ag. ''A^g ^Ag, ^Nr ^N^. ^N^ ^A^g. 1 1.35 
-A's -A„. -N^ -A's 2 1.28 


1 2 


: A^A^, AA, 2Ag ^Ag, ^Ns ^iVg, ^Ns ^As, ^iVg ^Ns 1 2.00 
^82(^02) : A^A, AA, ^Ns ^Ag, '^iVg '^A^g, '^iVg ^A^g 2 1.85 


■^P2 : iViV,iVA, AA, ^Ag ^Ag, ^iVg ^Ag, ^^Vg ^Ns, ^Ns ^iVg, 1 1.75 
^TVg ^Ag, ^TVg ^TVg 2 1.66 
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Fig. 2. The 5-wave phase shifts of A^A^ scattering. 

The model is also used to study the properties of 
deuteron. The results are given in Table 3. The bind- 
ing energy can be reproduced (we didn't fine time the 
strength of color confinement to get the better binding 
energy) as well as the D-w&ve component in the deuteron. 
However the root mean square radius is too small com- 
pared to experimental value. This shortcoming may be 



due to the fact that the deuteron wave function is only 
extending to 10 fm in this many channels coupling cal- 
culation. 
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Fig. 3. The mixing angles ci of A''A'' scattering. 
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Table 3. The properties of deuteron. 



recipes 


ac(MeV £m-2) 


EBiMoY) 


rms(fm) 


Pd 


1 


79.45 


1.0 


1.2 


4% 


2 


78.315 


2.2 


1.1 


4% 
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Fig. 4. The P-wave phase shifts of NN scattering. 
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Fig. 5. The D-wa,ve phase shifts of A'^A'^ scattering. 
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ThcTC' are problems remained and need to be studied 
further: (i) Up to now we don't know the confinement 
interaction between different color structures, this quark 
model study shows one need a phenomenological color 



confinement interaction different for different color par- 
tial waves but this might be artificial due to those approx- 
imations used in this calculation, (ii) The color screening 
effect due to unquench or qq excitation should be stud- 
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ied which might improve the above results. The effect 
of multi-body confinement interaction as obtained in the 
lattice QCD calculation should be studied too, because 
there are indications that it will effect the multi quark 
systems [llj. (iii) The n exchange is directly included 
in the Hamiltonian, is it possible to describe its effect 
directly from quark-gluon degree of freedom, as lattice 
QCD did ^6(1, by means of a more realistic nucleon inter- 



nal orbital wave function and a more sophisticated color 
confinement interaction as mentioned above. Such a cal- 
culation is quite involved numerically but it seems to be 
worth to devote. 
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